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Figure 2: Figure 2a 
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ABSTRACT 

We investigate the effects of heating generated by damping of Alfven waves 
on protostellar accretion disks. Two mechanisms of damping are investigated, 
nonlinear and turbulent, which were previously studied in stellar winds 
( |J atenco-Pereira fc Opher 1989a] , 1989b). For the nominal values studied, 
/ = 5v/va = 0.002 and F = w/fli = 0.1, where 5v, va and w are the amplitude, 
velocity and average frequency of the Alfven wave, respectively, and Q{ is the 
ion cyclotron frequency, we find that viscous heating is more important than 
Alfven heating for small radii. When the radius is greater than 0.5 AU, Alfvenic 
heating is more important than viscous heating. Thus, even for the relatively 
small value of / = 0.002, Alfvenic heating can be an important source of energy 
for ionizing protostellar disks, enabling angular momentum transport to occur 
by the Balbus-Hawley instability. 

Subject headings: accretion disks - MHD - stars: pre-main-sequence 
- waves 
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1. Introduction 

Accretion disks are a powerful energy reservoir and are expected to be found around 
pre-main sequence objects and active galactic nuclei and in binary systems. One of the 
most challenging problems that still remains concerning these objects is the nature of the 
mechanism by which they transport angular momentum. A very promising mechanism was 
proposed some time ago by Palbus fc Hawley (1991)1 (see also jBalbus fc Hawley 1998| ). 



based on a previous studied instability flVelikhov 1959| . |Chandrasckhar 196C ). An ambient 



magnetic field, together with a differentially rotating system, can give rise to a powerful 
instability that is very effective in transportating angular momentum in the disk (e.g., [Stone 
|et al. 1996|) . However, because of its magnetic nature, the instability needs a minimum 
ionization density, (n e /n n ) ^ 10~ 13 , in order to be effective (Gammie 1996). The presence 
of this minimum ionization density in protostellar accretion disks is uncertain. [Blaes fc| 
Balbus (1994)| made a linear study of the Balbus-Hawley instability (BHI) in a partially 



ionized medium. [Hawley fc Stone (1998)| did 3D numerical simulations and found that, in 



a disk that is not totally ionized, when 0.01 < ( / ypi/Q) < 100 (where jpi is the collision 
frequency and Q, the orbital frequency of the disk), there is angular momentum transport 
by BHI. Greater heating in the disk, in addition to that due to viscosity, will increase its 
degree of ionization and a larger fraction of the disk will be susceptible to the BHI. Our 
aim in this article is to analyze an extra source of plasma heating for the ionization of the 
neutral part of the disk, allowing the BHI to be effective in a greater part of it. 

A system that is magnetized and turbulent implies the generation of MHD waves. 
These MHD modes can be damped by several mechanisms, including collisions between 
ionized and neutral particles, resistivity, electron viscosity, electron thermal conductivity, 
ion viscosity ( [Kulsrud fc Pearce 1969| ) and mode coupling (Melrose 1980j , |Holzer et al. 1983] , 



Zwcibcl fc Josafatsson 1983). Damping of MHD modes can heat the gas. The damping of 



Alfven waves has been considered to be a heating mechanism for the solar corona (e.g., 
INarain fc Ulmschneider 1995|) and in the heating of filaments in cooling flows in galaxy 



clusters flFriag a et al. 199"7j ). It has also been considered as a source of acceleration and 



heating of the solar wind ( |Jatenco-Pereira fc Opher 1989a| , |Jatenco-Pereira et al. 19"94] ) 



as well as a generating mechanism of winds in late-type giant stars ( Holzer et al. 1983] , 
Jatenco-Pereira fc Opher 1989b|) . 



There may be other sources of energy which can heat the disk, in addition to viscous 
heating. For example, P'Alessio et al. ( 1999 )| showed that the irradiation from the central 
star, as well as ionization by energetic particles may be important for heating the disk. The 
irradiation is most important for R ^ 1 AU. As discussed by the authors, stellar irradiation 
heats the dust, while the viscous dissipation heats the gas. If the gas and the dust are well 
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mixed, the disk will have the same temperature for both components. Alfven waves will 
heat the gas, as will viscous dissipation. In principle, we can have these mechanisms (stellar 
irradiation, viscous dissipation and Alfvenic heating) all working together. 

In general, the structure of a protostellar accretion disk is thought to be as shown 
schematically in Figure |I[ Region 1 is the boundary layer between the disk and the star. 
This region can, in principle, extend up to the stellar surface. A difference in the rotational 
velocity of the two objects (disk and star) causes heating. The rotational velocity of the 
disk (~ 260 km s _1 ) is much greater than that of the star (~ 10 km s _1 ) and a large 
amount of energy is dissipated in a thin disk layer ( |Lynden-Bell fc Pringle 1974| ). However, 
if the star has a strong magnetic field (~ 1 kG; |Wang 19§Bj ), it can disrupt the disk, with 
the accretion proceeding along the field lines ( |Ghosh fc Lamb 1979a| , b; [Konigl 199 1| ). The 



magnetospheric accretion model accounts for many observational properties of T Tauri 
stars, such as the lack of a strong ultraviolet flux and the presence of blueshifted permitted 
lines ([Edwards 1997|) . A controversial point in this class of models is the precise location of 



the place where the magnetic field truncates the disk flHartmann 1998|) . Depending upon 



where this point (or region) is, the star can be spun up or down. Several wind models 
depend upon the determination of this location (e.g., |Shu et al. 1994j , 2000). Region 2 of the 
disk extends from ~ 0.01 AU to ~ 1 AU and is relatively close to the star. In this region, 
for 0.01 AU < R < 0.1 AU, the temperature is above 10 3 K, which is sufficient to couple the 
magnetic field to the gas, and it is likely that the BHI is present here. However, when the 
temperature falls below 10 3 K, the role of BHI is uncertain. Region 3 in the figure extends 
from 1.5 AU to ~ 5AU. |Gammie (1996)| claims that for R ^ 0.1 AU the disk separates 



into three regions. Two regions, above and below the midplane of the disk are ionized by 
cosmic rays ataining the necessary temperature to couple the magnetic field with the gas. 
Cosmic rays, however, do not penetrate deep enough to heat the midplane region and the 
magnetic field is thought to decouple from the gas preventing the BHI from operating. |lgca 



fc Glassgold (1999]] believe that X-rays from the central star are more effective in ionizing 



the disk than are cosmic rays. In region 4, the density of the disk is small enough for the 
disk to be totally ionized again. Somewhere in this region, the disk becomes optically thin 
and submillimeter observations provide information about the disk mass. 

In this article, the extra source of heating for ionizing the neutral part of a protostellar 
accretion disk is assumed to be due to the damping of Alfven waves. The damping 
mechanisms studied are nonlinear damping (section |2.1| ) and turbulent damping (section 
|2.2| ). For both damping mechanisms, we considered the cases of constant Alfvenic heating, 
exponentially varying density and temperature, as well as a layered model. We have made 
our calculations for the region that extends from 0.1 AU to ~ 100 AU. The results and 
discussions are presented in section [3], our conclusions are given in section f|. 
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2. Alfvenic Heating 



In order to quantify the amount of Alfvenic heating in a protostellar accretion disk we 
study the equations (in cylindrical coordinates) that describe a geometrically thin, optically 
thick, non-magnetized accretion disk flPringlc 1981 , Hartmann 1998 ): 



dF, 
dz 



9 GM, 

l vp ~W 



M 
3tt 



1- =^ 



R* 
R 



1/2 



(2) 



where F z is the energy flux in the vertical z direction, v is the viscosity, p is the fluid 
density, G is the gravitational constant, M* is the mass of the central star, R is the distance 
from the central star, £ is the surface density and R*, the radius of the central star. 



Integrating equation (II]), we obtain 



9 GM, 
-z/L — — . 
4 R 3 



(3) 



where we made use of the relation £ = J pdz. Substituting equation @ into equation @, 
we obtain the total flux of energy (erg cm~ 2 s _1 ) generated by viscous dissipation in the 
disk as a function of the disk radius: 



D 



Tot 



A 



3GM*M 



in 



1- — 



R> 
R 



1/2" 



(4) 



where M is the accretion rate of the disk. Equation @ is the gravitational potential energy 
liberated by the accretion flux. Since, in the present paper, we are examining the damping 
of Alfven waves as an extra source of energy for heating the disk, we have 



Drot — Dvis + Da, 

where Da is the energy generated by the damping of Alfven waves in the disk. 



(5) 



It is known that the nonlinear growth of the Balbus-Hawley instability results in 
turbulence ( [Stone et al. 1996|) and so, we can have MHD modes traveling throughout the 
disk. These MHD modes can be damped by several mechanisms (e.g., Holzer et al. 198 "3 



|Zweibel fc Josafatsson 1983|) . Here, we deal with two of these mechanisms: nonlinear 
damping and turbulent damping. These two damping mechanisms have been previously 



- 5 - 



investigated in a series of articles related to stellar winds by Jatenco-Pereira & Opher 
(1989a, 1989b) and Jatenco-Pereira et al. (1994). Here we investigate these damping 
mechanisms in relation to protostellar accretion disks. 

High amplitude waves cause nonlinear mode couplings and, in this process, the resulting 
modes are rapidly damped. In general, the waves that result from the damping processes 
are sound waves; thus, the nonlinear Alfvenic energy is, in general, converted into thermal 
energy. If the Alfven waves are not in phase, their energy is transfered both to the fast and 
slow magnetosonic modes (e.g., |Holzer et al. 198"3| ). In a high (3 plasma {(3 = 8irP/B 2 > 1, 



where P is the gas pressure and B 2 /8n, the magnetic field pressure), the mode coupling of 
the Alfvenic mode to the slow magnetosonic mode can be very strong because the phase 
velocities of the modes are essentially the same ([Melrose 1980[). 



Nonlinear damping of Alfven waves requires waves of great amplitude SB / B (where 
SB is the amplitude of the magnetic field of the wave and B is the strength of the magnetic 
field of the medium). The turbulent damping mechanism requires the presence of turbulent 
cells or vortices. If we have a turbulent medium in which SB/ B is large, in principle, we can 
have both mechanisms working together. In this paper we treat each mechanism separately. 

In the next sub-section, we consider the nonlinear damping mechanism of Alfven waves. 
(The damping rates used are obtained from Gongalves et al. (1998), hereafter GJO). 



2.1. Nonlinear damping of Alfven waves 

From equation (|5|), the total amount of heating in a magnetized accretion disk, assumed 
to be due only to Alfven damping and viscosity, is 

f+h r-+h 

D To t = D vis + / H A dz = D vis + 2 / H A dz, (6) 
J-h Jo 

where h is the half thickness of the disk, Ha = & w /La is the rate of Alfvenic heating, $„, 
is the Alfven wave flux, La = va/^w is the damping length, va = B/^Airp is the Alfven 
velocity and T w , the damping rate of the waves. The nonlinear damping rate is given by 

r 1 Wt„ f c s\ P(^ 2 ) m 

where £ is a constant with a value between 5 and 10, w is the characteristic frequency of 
the waves, c s is the sound velocity and p(Sv 2 ) j{B 2 /8ir), the energy density of the waves 
divided by the magnetic energy density ( |GJ(J| ). The energy density of the waves can also 
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be written as flWhang 19971) 



e A if v ,n = }:p(5v 2 ) + ^(5B 2 ). (8) 



The two terms on the right side of the equation @ are generally equal, so that 

tAlfven = P(8V 2 ). (9) 

Using the equations above, the wave flux can be expressed as 

$ w = p(Sv 2 )v A , (10) 
where (5v 2 ) is proportional to the degree of turbulence of the system. 

^From equations ((7|) and ([TOf) , we obtain the nonlinear Alfvenic heating rate: 

h nl = ^ 2 Wi / M 



The adiabatic sound velocity is c s = JjP/p and, for an ideal gas, P = p$tT/fi, where 
9ft is the universal gas constant. We define the parameter / by 

(5v 2 ) = fv\. (12) 
The average frequency w can be written as 

eB 

zu = FQ i = F , (13) 

rriiC 

where F(< 1) is a free parameter, f2j is the ion-cyclotron frequency and e, m ; and c are the 
electron charge, the ion mass and the velocity of light, respectively. 

Substituting equations ( JT^ ) and ( |TB"D in equation (|TTD, we have 

1/2 



8 mjC \ p J 
We then have, from equations @ and (0) 



(14) 
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D NL = 2 I" H NL dz = ^f'F^L V2 f\ p T)^BHz. (15) 

jo 4 m,c \ ll Jo 



rrtiC \ ji ) 



Gammie (1996) proposed a model in which a protostellar accretion disk has different 
layers with different ionization rates. Near the star (R <; 0.1AU), the temperature is greater 
than 10 3 K and, thus, collisional ionization is sufficient to couple the gas to the magnetic 
field. Beyond R = 0.1AU, the temperature drops and only small layers above and below 
the central plane of the disk (called "active layers") are able to sustain the Balbus-Hawley 
instability by cosmic ray ionization. The central layer, called the "dead zone", is neutral, 
with insufficient ionization to enable coupling. The active layers have a constant surface 
density, which depends on the cosmic-ray ionization rate, a poorly known quantity. |Gammie 



|(1996)| adopted the interstellar value QSpitzer fc Tomasko 1968|) for this rate and obtained 
S a = 10 2 g cm -2 , for the surface density of the active layers. Thus, if we assume that 
beyond the radius of 0.1 AU the surface density of the active layers is constant, we obtain 
from equation ([T5]) 

Dnl = ^lf F ^L ( ^V /2 S y 2 l h T l ' 2 B 2 z- l ' 2 dz, (R > 0.1AU), (16) 

4 TTijC \ fl J Jo 

where the relation £ = pz was used and it is assumed that the thickness of the central 
neutral region is very small compared to the thickness of the active layers (Gammie 1996). 

In quantifying the Alfvenic heating of a protostellar accretion disk, we consider various 
cases. The magnetic field in all of these cases, is held constant in the vertical direction, as 
was done in various previous papers (e.g., |Stone fc Norman 1994| ; [Armitage 199S| ; [Wardlc 



1999| ). First, we assume that the density and the temperature are constant throughout the 



disk (section |2.1.1| ). Then, we vary the density but still assume an isothermal disk (section 
|2.1.2| ). In the next case, we vary the temperature as well as the density (section |2.1.3| ). 
Finally (section |2.1.4j ), we consider the layered model of |Gammie (1996)| . In the following 
subsections, we assume an average effective value for / and F in the disk. 



2.1.1. Constant Nonlinear Alfvenic Heating 



Considering the density, temperature and magnetic field to be constant throughout the 
disk, we integrate equation (|15| ) and obtain 
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D 



NLconst 



4 rriiC 



1/2 



(17) 



where we used 7 = 1, since the temperature is constant everywhere. The temperature 
T d , adopted above, is obtained by assuming that the disk radiates as a blackbody and 
that the only source of dissipation is viscous heating. From D V i S = aT d (where a is the 
Stefan-Boltzmann constant), we obtain 



rp J 



3GM,M 
8ncr 



R 



R- 



(18) 



2.1.2. Nonlinear Alfvenic Heating with Exponentially Varying density 
According to the model of a stationary, geometrically thin accretion disk (|Pringlc| 



1981|) , if the temperature in the vertical direction is held constant, the resultant density 
profile is proportional to exp(— z 2 ). However, in a turbulent disk, it is expected that the 
decrease in density in the vertical direction is less dramatic, so that we assumed that the 
density falls as exp(-z). Thus, considering a disk with an exponentially varying density, 
p = p c e~ z / h , equation ( |T5| ) is given by 

1/2 



D NLpvar ^fF^_ (p c T d ) 1/2 B 2 h, 
5 miC \ p J 



(19) 



where we again used 7 = 1. 



2.1.3. Nonlinear Alfvenic Heating with Exponentially Varying density and Temperature 

In order to obtain a consistent temperature profile in the vertical direction, it is 
necessary to solve the equations that describe the disk in the z direction (e.g., P'Alessio 
1996| ; |D 7 Alessio et al. 1998| , 1999). This calculation requires a detailed treatment of the 



transport of the radiation and the advection of material by the BHI, which is beyond the 
scope of this paper. For simplicity, we have assumed that the temperature of the disk in 
the vertical direction falls off exponentially. So, in the case of a disk with both a density 
and temperature variation (T = T c e~ z ^ h ) the Alfvenic heating is given by 
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'NLvar — 2 q 



fF— f— V 2 (p c T c ) 1/2 B 2 h. (20) 
rriiC \ n ) 



2.1.4- The Layered Model 

In this case, to quantify the amount of Alfvenic heating of the disk, we adopt the 
model of |Gammie (1996]| . We consider only the region beyond 0.1 AU, where it is generally 



assumed that the disk is divided into three regions with different ionization densities. 
Adopting T = T c , together with a constant magnetic field, and integrating equation (|1 
we obtain 



D NLlay = ^f 4 F^_ (l®) 1/2 (SqTc h)1/2B2 (21) 
2 rriiC \ fj, ) 



2.2. Turbulent Damping 

A damping length, Lt, can describe the absorption of turbulent Alfven waves over a 
broad range of Alfven periods. Similar to a Kolmogorov spectrum, a turbulent cascade 
transfers wave energy from large to small scale regions. The dissipation of energy through 
turbulence is governed by the small scale linear absorption mechanism. Noting the similarity 
of Pb oc /c -5 / 3 and Kolmogorov turbulence in ordinary fluids^], the volumetric heating 
rate associated with the cascade is H T = p(5v 2 ) 3 ^ 2 / L corr , where L corr is a measure of the 
transverse correlation length (|Hollweg 1986| ). In terms of the damping length, 



L corr v(5v 2 ) 1/2 , (22) 



where L corr oc B x l 2 ( [Hollweg 1986| ) and v is the fluid velocity. 



In the present work, we treat Lt as an independent damping process, although it may 
be connected with non-linear damping. 

Then, the Alfvenic heating rate is given by ( |GJO| ) 



1 Where Pb is the observed magnetic power spectra in the solar wind and k is the wavenumber (Hollweg 
19871 ) . 
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va 

Using equation (|T^), we obtain 



H T = pB^fv\ = ( -j= \ BWp-W (24) 



In this case, the damping heating rate does not depend on temperature. Similar to our 
analysis in section |2.1| , we consider three cases in order to obtain the total dissipation due 
to turbulent damping of Alfven waves: 1) a disk with a constant density (section [2.2.1|) ; 2) 
a disk with a variable density (section [2.2.2|) and 3) a layered disk (section |2.2.3|) . 



2.2.1. The Constant Turbulent Damping 

We now consider a disk in which the density and the magnetic field are constant 
throughout. Integrating equation (|4]), we obtain 

D Tconst = 2 ( J=) ( ^) ' . (25) 



'Air , 



2.2.2. Variable Density 

Assuming that the density is varying exponentially with height (p = p c e~ z / h ), we can 
integrate equation (|24]), obtaining 

DT - = H-M{vT h - <26) 



2.2.3. Layered Disk 

Integrating equation fl24"|), using the layered model of constant surface density 
(p = Tl a /z) (section 2.1.4[ ), we obtain the total amount of heating for the turbulent damping 
of Alfven waves, 
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D Tlay = -\-L=\ — _ . (27) 



This equation is valid for the region beyond 0.1 AU, where it is assumed that there are 
active layers, created by cosmic ray ionization. 



3. Results and Discussions 

In Table [I], we show the values of the parameters of the disk that are used to calculate 
the different types of Alfvenic heating for region 2, except for those related to the layered 
model of |Gammie ( 1996J] (D^Liay and D T i ay ). These values (except the magnetic field) are 



obtained using a model of a stationary, geometrically thin, optically thick, a accretion disk 
(e.g., |Pringle 1981| , [Frank, King fc Raine 1992] , [Shakura fc Sunyaev 1973]) with the opacity 



law of [Bell fc Lin (1994~| in the iron-silicate regimen. The magnetic field is obtained from 



the model of Gammie (1996 



In Table |2j, we show the disk parameters obtained from the model presented by 
|Gammie (1996)| for region 2. With the parameters shown in Table ||, we calculate the 



layered Alfvenic heating for region 2. The values of both models are obtained considering a 
disk surrounding a T Tauri star of 0.5Mq, which is accreting mass at a rate of 10 -8 Mq 
year -1 . We adopt a value of 0.6 for the mean molecular weight /x. 

Temperatures associated with each dissipated energy mechanism may be calculated, 
assuming that the disk radiates like a blackbody. Thus, there is the following temperature 
associated with D NLconst , for example: 

/ n \ V 4 

rp I ^NLconst \ (OQ\ 
1 const — I — I > {^O) 

where T const is the temperature and a, the Stefan-Boltzmann constant. 

In Figure 2, we plot the temperatures associated with nonlinear Alfvenic dissipation 
for disk region 2, calculated from equation (f28|) with / = 0.002 (eq. fl2]|), F = 0.1 (eq. 
flcj ) and £=5. Figure 2a shows the temperatures associated with each kind of dissipation. 
T vis is obtained from the viscous dissipation of the stationary, geometrically thin accretion 



In this regime, the opacity is dominated by iron and silicate grains. The corresponding range of 
temperature is between ~ 10 3 K to 203 K (Gammie 1996). 
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disk model (Table [T]). For comparison, we plot T irr , the temperature associated with the 
model of the irradiated disk ( P'Alessio et al. 1998| ). Plotted values in the figure were 
obtained from Figure 3 of D'Alessio et al. (1998). The temperature calculated using 
Alfvenic dissipation in the layered model (Ti ay , shown as circles in Fig. 2a) is highest for 
R > 0.2 AU. The other temperatures (T const , T var , T pvar ) do not differ from one another 
very much. T vis = (D^/cr) 1 / 4 decreases with increasing radius, becoming less than the 
other temperatures for R > 0.4 AU. The temperature obtained from the irradiated model 
is greater than T const , T pvar , T var for all radii considered. We plot the total temperatures 
obtained from our models together with T vis in Figure 2b. They are calculated from the 
following relations: 



T totl = [T d 4 + T t 



consti 



,1/4 



rp _ \rpA I rpA ll/4 _ 

totl — I 1 d + 1 pvarl ~ 
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(29) 
(30) 
(31) 
(32) 



We note the increase in temperature. Highest temperatures are again obtained when 
the layered model is used. The different density and temperature profiles, used to obtain 
the other Alfvenic dissipations {D^Lconst, D^ipvar and D^L var ), do not affect the final total 
temperature obtained. They are practically equal for all radii. 

In Figure 2c, we plot the total temperatures, taking into account the temperature from 
the irradiated model. The increase in temperature is much more significant than when only 
viscous dissipation is considered. Again, the density and temperature profiles do not alter 
the final temperature apreciably. The highest temperature is obtained from the layered 
model. T irr dominates the total temperature values, T T1 , T T2 , T T3 , which indicates that 
irradiation may be more important than the Alfvenic heating mechanisms which produced 
T C onst, Tpvar, T var for this part of the disk. 

We further extend our calculations to include a larger amount of the protostellar disk. 
We divide the disk according to the different opacities laws, as given by |Gammie (1996 



We make our calculations for two additional regions (regions 3 and 4 of Figure 1). region 
3 extending from 1.5 AU to 5.0 AU and region 4 extending from ~ 5 AU to ~ 100 AU. 
The opacity law for region 3 is given by k = 2 x 10 16 T~ 7 cm 2 g -1 and for region 4, by 
k = 2 x lO-^cmV 1 (Bell & Lin 1994). 
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In Figure 3a, we show the temperatures obtained from nonlinear heating models, 
calculated for region 3 (1.5 AU < R < 5 AU). For this region, T vis (showed again as up 
triangles in figure) is the smallest temperature. The temperature obtained for the layered 
model (circles) is greater than that obtained for other models for all radii considered. The 
photospheric temperature obtained from the irradiated model is compared with T var (stars 
and squares in the figure, respectively). In this figure, we note a small separation between 
the different profiles of temperature and density. 

In Figure 3b, we plot the temperatures obtained from nonlinear heating models for 
region 4, which extends from ~ 5 AU to ~ 100 AU. We note small changes compared with 
Figure 3a. We can observe an inversion of the order of the temperatures. For this case, T var 
is less than T const and T pvar . However, this inversion is not very dramatic. T[ ay , again, is 
the highest temperature and Tj rr does not differ very much from the Alfvenic temperatures 
T const and T pvar . In this part of the disk, the presence of layers with different degrees of 
ionization is less certain. Due to the small values of density, the penetration of cosmic rays 
can be deeper, perhaps, sufficiently to ionize this entire disk region. 

In Figure 4a, we show the turbulent Alfvenic heating and viscous dissipation. For all 
radii, Alfvenic heating is greater than the viscous heating. This means that the turbulent 
mechanism for damping can be much more efficient to dissipate energy than the viscous 
mechanism. Moreover, we can see that, in this case, the turbulent Alfvenic heating 
calculated, using the layered model D T i ay (eq. p7|), is less than both the Alfvenic heating 
with constant density Dxconst (eq. [[25[]) and the Alfvenic heating with exponentially variable 
density Dxvar (eq. [pE] ). However, DTconst and Dxvar are practically the same for all radii. 
This means that turbulent damping is not very dependent on density. In Figure 4b, we 
plot the turbulent Alfvenic heating and viscous dissipation for region 3 (1.5 AU < R < 
5.0 AU). We note that viscous dissipation is very small, when compared with turbulent 
dissipation. For D pvar and D const , there is a tendency for an increase at the beginning but 
after R ~ 3 AU this tendency is reverted and the energy decreases. Again, D[ ay is not very 
significant, when compared with the two other turbulent dissipations. In Figure 4c, we 
show the turbulent dissipated energies for region 4 (5 AU ^ R ^ 100 AU). Again, turbulent 
dissipation is greater than viscous dissipation. 

We know that an important parameter for disk models is the electron fraction, 
x = n e /n n , where n e is the number density of electrons and n n , the number density of 
the neutral specie. |Gammie (1996)| calculated an upper limit for this fraction, assuming 



ionization by cosmic rays. He obtained 
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x ={-^) ' =1 - 6xl0_12 (^) (TT^fcr) [t^—,) > (33) 

where £ is the cosmic-ray ionization rate, /3 is the recombination coefficient and uh, the 
number density of hydrogen. This expression is valid for T < 10 3 K. Using equation (|33| ) 
and Td, E and h from Table 0, we find that the electron fraction increases with radius. 
However, in order to evaluate x, [Stone et al. (200UJ| used an expression derived from the 



Saha equation, which took into account the thermal ionization of potassium: 

/oiy in 15 \ 1/2 

x = a 1 ^ 4 £±ii_^_ exp (_50, 370/2T), (34) 
\ n n ) 

where a is the abundance of potassium. In this case, x decreases with radius (The expression 
is valid for T < 2000/^). It is possible that, for the range of radii and temperature considered 
in this paper, the expression used by pammie (1996)| for x is the most likely to be valid. 



However, the determination of x is still very uncertain. 



4. Conclusion 

In this paper, we analysed the important role that the damping of Alfven waves 
may possibly play in the inner and intermediate regions of protostellar accretion disks 
(0.1 < R(AXJ) <z 100). We considered two damping mechanisms: nonlinear and turbulent. 

We considered various profiles of density and temperature, with a constant magnetic 
field, in quantifying the heating generated by the damping of the waves. Concerning 
nonlinear damping, in the first case, the density and temperature are held constant 
throughout the disk. The resultant heating is called D^iconst- In the second case, the 
density is exponentially varying, while the disk is isothermal (DNipvar)- Both temperature 
and density are exponentially varying in the third case, Dnlvht- The heating D^Liay, 
the fourth case, is obtained when nonlinear damping of Alfven waves occurs in the disk 
model of |Cammie (1996")| (the layered model). We note that D NLlay is greater than the 
viscous dissipation, D vis for R > 0.2AU. We see that, for R > 0.5 AU all the forms of 
nonlinear Alfvenic heating are greater than viscous heating. When the damping occurs by 
turbulent mechanisms, all the resultant turbulent Alfvenic dissipations are greater than 
viscous dissipation for all considered radii. Figure || shows the calculated temperatures for 
our model for region 1 (0.1 AU < R < 1.4 AU). We note an increase in temperature for 
this region in Figure 2b. Thus, we conclude that, when viscous dissipation is insufficient to 
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ensure the necessary degree of ionization for BHI to occur, the damping of Alfven waves can 
be an alternative source of energy for ionization. When taking into account the irradiation 
from the central star, the increase in temperature is much greater than when only the 
viscous dissipation is considered. Alfvenic heating in the layered model is very significant, 
greater than all other dissipations considered. 

In Figures 3a and 3b we show the temperatures obtained from the nonlinear heatings 
for the outer parts of the disk. Again, the same order for the temperatures remains, with 
little differences shown. The viscous dissipation becomes smaller and less significant for 
increasing radii. The irradiation is a rather important source of heating, comparable with 
nonlinear Alfvenic D const , D pvar and D var . In Figures 4b and 4c we see that turbulent 
dissipation is very efficient for heating the disk. We thus conclude that Alfvenic heating 
can be important in the outer parts of the disk. Considering only viscous dissipation, the 
temperatures reached are not sufficient to ensure the development of BHI and, consequently, 
of turbulence. If we take into account the irradiation from the central star, temperatures 
can reach values which ensure the necessary coupling between the gas and the magnetic 
field. 

In the present paper, it is to be noted that we used /(= Sv/va) = 0.002 and 
F(= w/Qi) = 0.1. Nonlinear Alfven heating is proportional to f 4 F and turbulent Alfven 
heating, to f 3 . Increasing (decreasing) / and F from these nominal values increases 
(decreases) the Alfven heating. We have shown, in this paper, that even for the relatively 
small value of / = 0.002, Alfvenic heating can be important in protostellar disks. 

M.J.V. would like to thank the Brazilian agency FAPESP for financial support. V.J.P 
and R.O. thank the Brazilian agency CNPq for partial support. The authors would like to 
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Fig. 1. — Sketch of a protostellar accretion disk. Region 1 corresponds to the boundary 
between the disk and the star. Region 2 extends from ~ 0.01 AU to ~ 1 AU. Region 3 
extends from ~ 1 AU to ~ 5 AU. Beyond this region, region 4, the disk becomes optically 
thin. We make our calculations for regions 2,3, and 4, extending from 0.1 AU to ~ 100 AU. 
(See text for a more detailed discussion.) 

Fig. 2. — Temperatures associated with nonlinear Alfvenic dissipation, calculated for 0.1 
AU < R < 1.4 AU, assuming that the disk radiates like a blackbody. 2a) Each curve in 
this figure describes the temperature dependence vs radius calculated using the dissipation 
energy mechanism indicates (eq. ^8|).The values attributed to T irr , the temperature for an 
irradiated disk model was extracted from Figure 3 of D'Alessio et al. (1998); 2b) In this 
figure, we show the total temperatures reached in our models, taking into account only the 
viscous dissipation. They are obtained by summing the energy radiated by corresponding 
to viscous dissipation with the Alfvenic heatings (equations 29-32). 2c) This figure shows 
the total temperatures summing viscous dissipation, the irradiated model and each Alfvenic 
heating. 

Fig. 3. — Temperatures associated with nonlinear Alfvenic dissipation, calculated for: 3a) 
for region 3 (1.5 AU < R < 5 AU) and 3b) for region 4 (5 AU < R ~ 100 AU). 

Fig. 4. — Dissipated energies (in units of erg cm -2 s _1 ) as a function of disk radius (in AU), 
obtained from the turbulent heating model, for: 4a) region 2 (0.1 < R < 1.4 AU), 4b) region 
3 (1.5 AU < R < 5 AU) and 4c) region 4 (5 AU < R ~ 100 AU). We observe that turbulent 
Alfvenic heating is much greater than viscous heating for all density profiles. 
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Table 1. Parameters of a stationary, geometrically thin, optically thick protostellar 
accretion disk, with the opacity law of Bell & Lin (1994) in the iron-silicate regime. 



R (AU) 


Tc (K) 


T d (K) 


h (AU) 


£ (g cm- 2 ) 


p c (g cm 3 ) 


0.1 


1362.36 


369.28 


0.008 


66.92 


5.31 x 10~ 10 


0.2 


718.80 


226.31 


0.017 


50.61 


1.95 x 10~ 10 


0.3 


489.80 


169.04 


0.026 


42.47 


1.08 x 10~ 10 


0.4 


371.99 


137.20 


0.035 


37.36 


7.09 x 10~ 9 


0.5 


300.09 


116.60 


0.044 


33.77 


5.11 x 10" 9 


0.6 


251.60 


102.05 


0.053 


31.06 


3.90 x 10~ 9 


0.7 


216.67 


91.15 


0.062 


28.92 


3.11 x 10~ 9 


0.8 


190.29 


82.63 


0.071 


27.18 


2.55 x 10~ 9 


0.9 


169.67 


75.78 


0.080 


25.73 


2.14 x 10~ 9 


1.0 


153.09 


70.12 


0.089 


24.48 


1.83 x 10~ 9 


1.1 


139.48 


65.36 


0.098 


23.41 


1.59 x 10~ 9 


1.2 


128.10 


61.30 


0.107 


22.47 


1.40 x 10~ 9 


1.3 


118.45 


57.78 


0.116 


21.63 


1.24 x 10~ 9 


1.4 


110.15 


54.71 


0.126 


20.88 


1.11 x 10~ 9 



Note. — Here, R is the disk radius, T c is the temperature of the midplane of the disk, Td is the temperature 
due to viscous dissipation, h is the disk semi-thickness, B is the magnetic field (obtained from the model of 
Gammie 1996) and p c , the midplane density. 
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Table 2. Parameters of the disk calculated using the layered model of Gammie (1996), for 

a region extending from 0.1 AU to 1.4 AU. 



R (AU) 


T c (K) 




h (AU) 


B(G) 


Pc (g cm 3 ) 


0.1 


1005.06 


302.64 


0.003 


4.96 


3.51 x 10~ 9 


0.2 


663.09 


210.32 


0.008 


2.66 


1.53 x 10~ 9 


0.3 


519.90 


169.99 


0.012 


1.84 


9.40 x 10~ 10 


0.4 


437.48 


146.16 


0.017 


1.42 


6.66 x 10~ 10 


0.5 


382.66 


130.01 


0.023 


1.16 


5.09 x 10~ 10 


0.6 


343.01 


118.14 


0.028 


0.99 


4.09 x 10~ 10 


0.7 


312.70 


108.95 


0.034 


0.86 


3.40 x 10~ 10 


0.8 


288.63 


101.58 


0.040 


0.76 


2.90 x 10~ 10 


0.9 


268.93 


95.49 


0.046 


0.69 


2.51 x 10~ 10 


1.0 


252.46 


90.35 


0.052 


0.62 


2.22 x 10~ 10 


1.1 


238.43 


85.94 


0.059 


0.57 


1.98 x 10~ 10 


1.2 


226.30 


82.10 


0.065 


0.53 


1.78 x 10~ 10 


1.3 


215.69 


78.72 


0.072 


0.49 


1.62 x 10~ 10 


1.4 


206.31 


75.72 


0.078 


0.46 


1.48 x 10~ 10 



Note. — In this tabic, R is the disk radius, T c is the temperature of the midplanc of the disk, Td is the 
temperature due to viscosity dissipation, h is the semi-thickness of the disk, B is the magnetic field and p c , 
the midplane density. 



